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Quantum Markov networks are a generalization of quantum Markov chains to arbitrary 
graphs. They provide a powerful classification of correlations in quantum many-body systems — 
complementing the area law at finite temperature — and are therefore useful to understand the pow- 
ers and limitations of certain classes of simulation algorithms. Here, we extend the characterization 
of quantum Markov networks [151 123] and in particular prove the equivalence of positive quantum 
Markov networks and Gibbs states of Hamiltonians that are the sum of local commuting terms on 
graphs containing no triangles. For more general graphs we demonstrate the equivalence between 
quantum Markov networks and Gibbs states of a class of Hamiltonians of intermediate complexity 
between commuting and general local Hamiltonians. 



I. INTRODUCTION 

A Markov network is an object used to represent the 
correlation structure of a probability distribution over a 
very large number of random variables. In this repre- 
sentation, random variables are located at the vertices 
of a graph and the edges encode correlations between 
neighboring vertices. While distant vertices can be cor- 
related, their correlations are mediated by intermediate 
random variables located along the paths connecting the 
two vertices; distant variables are conditionally indepen- 
dent given the variables separating them. Since a se- 
quence of conditionally independent variables defines a 
Markov chain, Markov networks are their natural gener- 
alization from chains to arbitrary graphs [TU [TH1 [HI H2] • 

An important motivation for the characterization of 
Markov networks stems from the existence of powerful 
algorithms and heuristics to solve or approximate the so- 
lution to statistical inference problems on these graphi- 
cal models, see, e.g., H3 [55] and references therein. 
The archetypical inference problem consists in determin- 
ing the marginal probability of a given random variable 
after having observed some of the other random variables. 
Applications of such inference problems range from med- 
ical diagnostics [25] to digital communication [20], and 
from computer vision |17) to gene expression analysis [5] . 

Another application, which will be the main focus in 
this article, is statistical physics. In this setting, parti- 
cles (e.g. Ising spins) are subjected to local interactions 
on a regular lattice or a more general graph, and we are 
interested in computing some correlation function at a 
given temperature. The thermal equilibrium state of the 
system is a Gibbs probability distribution, which is a 
function of the interaction between particles. A pow- 
erful characterization of such graphical models, due to 
Hammersley and Clifford [5] , shows that the set of Gibbs 
distributions arising from local interactions on a given 
graph coincides with the set of Markov networks on that 
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FIG. 1: The Markov condition demands that the mutual in- 
formation between region A and the rest of the system C, 
condition on its boundary B, be zero. More formally, this 
condition I(A : C\B) — must be satisfied whenever A, B, 
and C are disjoint regions such that all paths from A to C go 
through B. 



graph. 

In this article, we study quantum mechanical version of 
this characterization. Quantum Markov networks, intro- 
duced in |15) . can be defined similarly to classical ones, 
in terms of conditional independence. Informally, we say 
that a many-body quantum state is a Markov network 
when the mutual information between a region A and 
the rest of the system, conditioned on the boundary of 
region A, vanishes, see Fig. [T] Our general goal is to un- 
derstand the relation between such states and quantum 
Gibbs states. The non-commutativity of quantum oper- 
ators is the main obstacle to a direct generalization of 
the classical Hammersley-Clifford Theorem. 

It was previously known [15] that quantum Markov 
networks can be expressed as Gibbs states of local Hamil- 
tonians, but that the converse is not generally true. 
While the converse holds whenever the terms in the 
Hamiltonian mutually commute — i.e., local commuting 
Hamiltonians, as we call them for short, give rise to 
Gibbs states that are Markov networks — it was unknown 
whether this commutation condition is necessary. That 
is, whether there is an equivalence between positive quan- 
tum Markov networks and Gibbs states of commuting 
Hamiltonians. We show, as our main result, (Theorem|4|, 
that there exists such an equivalence between local com- 
muting Hamiltonians and quantum Markov networks on 
graphs that contain no triangular cells, extending a re- 
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suit previously derived for cycle- free graphs [53] . We also 
construct an example of a positive quantum Markov net- 
work on a triangular lattice, which is not a Gibbs state 
of a commuting Hamiltonian. 

The motivation for our work stems from several 
sources. Quantum many-body systems are believed to 
obey an entanglement area law at zero temperature: in 
the ground state of a gapped interacting quantum sys- 
tem, the entanglement between a region A and the rest 
of the system should only increase proportionally to the 
size of the boundary of A, not its volume. The area law 
reflects the fact that the quantum correlations in the sys- 
tem are short ranged. At finite temperature, it has been 
shown |26] that the mutual information between region 
A and the rest of the system obeys an area law. This 
finite-temperature characterization does not describe en- 
tanglement, however, since both classical and quantum 
correlations are picked up by mutual information. Con- 
ditional mutual information could be an important tool 
to study entanglement area laws at finite temperature 
because its vanishing demonstrates that quantum corre- 
lations are short range, while classical correlations are 
Markovian but can extend to arbitrarily long range. 

Thus, the Markov condition ensures that the quantum 
many-body state has only finite range entanglement. In 
a system with topological order, the conditional mutual 
information would be zero for any topologically trivial 
region A, however, there must exist a non-trivial region, 
for instance when A is a ribbon wrapping around a torus, 
such that the mutual information between A and the rest 
of the system, conditioned on the boundary of A, is equal 
to the topological entanglement entropy of the system 
[TU1 HI)] . The state obtained from the uniform mixture of 
all ground states of a topologically ordered system, such 
as Kitaev's toric code [12], can form a Markov network 
however. 

Another goal of this line of research is to leverage the 
power of the various algorithms and heuristics that have 
been developed in the classical setting to the quantum 
realm. For instance, belief propagation can be used to 
solve statistical inference problems involving many cor- 
related random variables. It is exact when the underly- 
ing statistical model is a Markov network on a cycle-free 
graph, and works remarkably well in other contexts. Be- 
lief propagation and its variational dual have been gen- 
eralized to the quantum setting [7J (T3] [TSl [53] ■ While 
numerical experiments have shown that they can provide 
reliable approximations, the present work could help un- 
derstanding their validity more rigorously. 

Lastly, our work is motivated by questions in quan- 
tum Hamiltonian complexity. One central question in 
this field of research is to determine the computational 
complexity of estimating the ground state energy of a 
quantum many-body system, or equivalently its free en- 
ergy at finite temperature, for which belief propagation 
can be used. A special case of this question concerns 
local commuting Hamiltonians. Although at first glance 
commuting Hamiltonians seem to be entirely classical, 



they can exhibit important quantum phenomena such as 
topological order [12]. It remains unclear at this stage 
whether the ground state of commuting Hamiltonians 
embody the same complexity as the ground state of arbi- 
trary local Hamiltonians. This problem is the subject of 
intensive studies, see e.g. [4j [6l [24] , and some of the ques- 
tions addressed in this article are directly related to it. 
In particular, quantum Markov networks seem to offer an 
intermediate structure between commuting and general 
local Hamiltonians. 



II. QUANTUM MARKOV NETWORKS 
A. Quantum Markov chains 

Before defining quantum Markov networks, we review 
the more familiar concept of quantum Markov chains. 
Throughout this article, we consider a finite dimensional 
system composed of TV particles labeled a = 1, 2, . . . , N, 
each associated to a finite-dimensional Hilbert space Tl a . 
The Hilbert space of the entire system is H = ^ a T~L a - 
The state of the system is a non-negative matrix p : % — > 
% with trace equal to 1. 

Conditional independence will play a crucial role in the 
definition of a Markov network. Given a density matrix 
Pabc on 3 disjoint collections of particles labeled A, B, 
and C, the conditional mutual information between A 
and C given B is defined as 

I(A : C\B) = S(AB) + S(BC) - S(ABC) - S{B) (1) 

where S(p) = — Trplogp is the von Neumann entropy, 
and S(X) = S(px) denotes the von Neumann entropy 
associated to the density matrix px = Tr^o of subsystem 
X, obtained by taking the partial trace over the comple- 
ment X of X. We say that A and C are independent 
conditioned on B if and only if I {A : C\B) = 0. 

Conditional independence implies that A-B-C forms 
a quantum Markov chain, i.e. it is possible to create 
the state Pabc starting with the state pab ® |0)(0|c and 
applying a physical transformation that involves the sys- 
tems B and C only [5]. 

A useful characterization of independence of A and C, 
conditioned on B [8] , is that the Hilbert space of subsys- 
tem B decomposes as, 

H B = Y. U -'-B^A®Hj:B^C (2) 
J 

such that pabc = J2jPJ a J-AB ® oj-.bc where pj is a 
probability distribution and each uj-ab and <jj-.bc & re 
density matrices with support only on Hj-b^-a ® 7~La 
and Hj-.b^c ®Hc respectively. 

We find it convenient to express the above characteri- 
zation as the following decomposition of the state into a 
product of commuting operators, 

Pabc = AabAsc with [A A b , A B c] = (3) 
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and the A's are non-negative. In this equation, we make 
use of a notation where Qx is an operator that acts only 
on the Hilbert space of X and is trivial everywhere else. 
It should be understood that Qx is an operator on X 
tensor product with the identity on the complement of 
X, but we omit these identity matrices throughout to 
simplify the notation. When the density matrix is posi- 
tive pabc > 0, we can take the logarithm of Eq. ^ and 
arrive at 



Pabc = e HAB+HBC 



ith [H AB ,H BC ] -0, (4) 



which is the standard expression for the Gibbs state 
associated to the Hamiltonian H = Hab + Pbc ( we 
rescale the Hamiltonian by —1 /temperature throughout 
to lighten the notation). 

The above construction generalizes to more than 3 sys- 
tems. Consider a collection of N quantum systems with 
density matrix p, and suppose that 



1(1,.. .k- 1 : k + l,...,N\k) = 



(5) 



for all k = 2, 3, . . . , N — 1. These conditions can be un- 
derstood by imagining that the particles are arranged in 
a chain, and that conditioned on any site k, the sites to 
the left of k are independent of the sites to the right of k. 
For each site, p decomposes according to Eq. Q. As will 
be shown in Sec III C| the decompositions for subsequent 
sites may be iterated, arriving at the conclusion that 
p = e H where H = J2k hk,k+i with [h k<k+1 , h k >,k'+i] = 0. 
Thus, quantum Markov chains are Gibbs states of local 



commuting Hamiltonians. As we will see in Sec. Ill B , the 
converse implication also holds and we have a complete 
equivalence between positive quantum Markov chains 
and local commuting Hamiltonians in ID. 



B. Quantum Markov networks 

A quantum Markov network is a generalization of 
the above construction from a chain to an arbitrary 
graph. A graph G is composed of a set of vertices V 
and edges E. The particles are located on the vertices 
V = {1,2, ... , N} of a graph G = (V,E); we use the 
same label a for a vertex and the particles located on 
it. Let A, B, and C be three disjoint subsets of V . We 
say that B shields A from C whenever all paths on G 
beginning at a vertex in A and ending at a vertex in C 
pass through a vertex in B. For instance, the boundary 
of a region A — composed of all the vertices that share an 
edge with a vertex in A but do not belong to A — always 
shields A from the rest of the lattice, as illustrated on 
Pig-H 

We say that the pair (p, G) , formed by a graph G with 
N vertices and a density matrix p for the N particles lo- 
cated on the graph G, is a quantum Markov network when 
I (A : C\B) = for all disjoint subsets A, B, C C V such 
that B shields A from C. In other words, in a Markov 
network, A-B-C forms a Markov chain when all paths 



from A to C go through B. Classical Markov networks 
are defined similarly, but with the Shannon entropy re- 
placing the von Neumann entropy and a probability dis- 
tribution replacing the density matrix. 

In Ref. |15j . it was shown that positive quantum 
Markov networks — those with p > — are Gibbs states 
of local Hamiltonians. For a general graph G, a local 
Hamiltonian is defined to be of the form 



H 



Qec 



(6) 



where C denotes the set of cliques of the graph G. Recall 
that a clique of G is a complete (fully-connected) sub- 
graph of G. For instance, on a regular square lattice, the 
Hamiltonian would be a sum of terms acting on nearest 
neighbors, while on a triangular lattice it would also in- 
clude three-body terms acting on the sites of a triangular 
cell. 

This characterization leaves open the converse impli- 
cation, namely under what circumstances do quantum 
Gibbs sates form Markov networks. This question is the 
focus of the following section. 



III. QUANTUM GENERALIZATIONS OF THE 
HAMMERSLEY-CLIFFORD THEOREM 

We begin by formally stating the classical characteri- 
zation Theorem: 

Theorem 1 (Hammersley and Clifford |5j) Let 

G = (V, E) be a graph and P(V) be a positive probability 
distribution over random variables located at the vertices 
of G. The pair (P(V),G) is a positive Markov network 
if and only if the probability P can be expressed as 
P(V) 



H{V) = h Q (Q) 
Qec 



(7) 



is the sum of real functions Hq(Q) of the random vari- 
ables in cliques Q, and Z is a normalization constant. 

This Theorem establishes a complete equivalence be- 
tween positive classical Markov networks and Gibbs dis- 
tributions arising from local Hamiltoinians, and in the 
next sections we will partially generalize this equivalence 
to the quantum setting. 



A. Markov implies locality 

One direction of this Theorem holds in the quantum 
setting: 

Theorem 2 (Leifer and Poulin [15]) Let G = (V,E) 
be a graph and p be a density matrix for the particles 
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located at the vertices of G. If the pair (p,G) is a pos- 
itive quantum Markov network, then the state p can be 
expressed as p = e H where 

H=J2 h Q ( 8 ) 

Qec 

is the sum of Hermitian operators hg on the particles 
located in cliques Q. 

Here, we present a proof that is much simpler than the 
one presented in [TS], although similar in essence. The 
key idea is to make a cumulant expansion of the "effective 
Hamiltonian" , formally defined as H — log p. Recall that 
any operator H acting on N particles can be uniquely 
decomposed into a sum 

H = J2 K * ( 9 ) 

xev 

where X runs over all the subsets of N particles. What 
makes this decomposition unique is that each cumulant 
Kx obeys TryKx — for any Y C X, i.e. they are 
"partial-traceless" . This property also implies that the 
cumulants are orthogonal with respect to the Hilbert- 
Schmidt inner product Tt(KxKy) = SxyTrKjci an( i it 
follows that Tr(HKx) = TrK x . 

Given these properties, we can easily prove the theo- 
rem. All that is required is to show that Kx = when- 
ever X is not a clique. Consider such a region X. Then, 
there exists two vertices a and c £ X that are not linked 
by an edge, and so the region B = V — a — c shields a 
from c. Repeating the arguments leading to Eqs. ([3j|4]), 
we arrive at H = H a s + Hb c , so 

Tt(HK x ) = Tr(H aB K x ) + Tr(H Bc K x ) 

= Tr(H aB [Tr c K x ]) + Tr(ff Bc [Tr a X x ]) = 

by the partial-tracelessness property of cumulants. This 
implies Trif^ — 0, and so Kx — as claimed, which 
completes the proof. 

B. Locality and commutativity implies Markov 

Unlike the classical case, it is not true, however, that 
an arbitrary Gibbs state of a local Hamiltonian Eq. ^ 
yields a positive quantum Markov network, as one can 
easily find examples of local Hamiltonians that do not 
generate Markov networks, even in ID |15j . Specifically, 
any Hamiltonian that cannot always be decomposed into 
a sum of commuting terms Hab + Hbc across region 
B that shields A from C generates a Gibbs state which 
violates condition Eq. Q. On the other hand, when all 
of the terms in the Hamiltonian commute, Eq. Q is al- 
ways satisfied, implying that the Gibbs state of any local 
commuting Hamiltonian is a quantum Markov network. 

Theorem 3 Let G = (E, V) be a graph and H = 
Sqgc ^Q> IAqj^Q'] = 0, be a local commuting Hamil- 
tonian on that graph. Then (p, G) is a positive quantum 



Markov network for p = -^e H , where Z is a normaliza- 
tion constant. 

To prove this theorem, we can first consider a restricted 
set of partitions A, B, and C — V — A — B such that B 
shields A from C. In this case, we can easily show that A- 
B-C forms a Markov chain from the observation that no 
clique Q can overlap simultaneously with region A and C, 
otherwise they would not be shielded by B. Thus, we can 
(non-uniquely) assign each clique to either region AB or 

BC and decompose H = J^QeAB + Sqsbc ^Qi ^ wo 
terms that obviously commute, ensuring that the Gibbs 
state is a Markov Network, c.f. Eq. Q. 

We can readily extend the result to an arbitrary choice 
of regions A and C that arc shielded by B, by the follow- 
ing procedure: expand regions A and C to AA' and CC' 
until AA'BCC' fills the entire lattice, and B still shields 
A A' from CC . This can be done by recursively expand- 
ing region A to include all of its neighbors that are not 
in B, and similarly for C. At the end of this process, 
any isolated islands that are not included in either the 
expanded A, the expanded C, or B can be included in A' . 
Then, because B shields A A' from CC = V-AA'-B, it 
follows from the argument of the previous paragraph that 
AA'-B-CC is a Markov chain, so I (AA' : CC'\B) = 0. 
Then— using I(AA' : CC'\B) > I (A : C\B), which is 
a straightforward consequence of strong sub-additivity 
of entropy [18] — we obtain that I(A : C\B) — 0, so 
A — B — C is also a Markov chain. 

We also note that the above result may be extended 
to non-positive density matrices that are projectors onto 
eigenspaces of arbitrary commuting Hamiltonians, in- 
cluding for example the projector onto the code space of 
any local stabilizer code. Since any such density matrix 
may be written as a product of the projectors onto the 
corresponding eigenspace of each term in the Hamilto- 
nian, it follows from Eq. ^ that conditional indepen- 
dence is satisfied for any set of partitions A, B, and 
C = V - A - B such that B shields A from C. Con- 
ditional independence for arbitrary choice of regions A 
and C that are shielded by B, then follows similarly from 
strong sub-additivity of entropy, implying that any such 
density matrices are quantum Markov networks. 

C. When does the Markov condition imply 
commutat ivity ? 

As seen in the proof of Theorem [3j a quantum Markov 
network can be defined by demanding that the condition 
I (A : C\B) = holds for any region A, B and C that 
span all vertices with B shielding A from C; the condition 
for regions A, B, and C that span only a subset of the 
vertices follow from strong subadditivity. Equation Q 
provides a general characterization of such states, which 
motivates the following definition. 

Definition 1 (Shield commuting Hamiltonian) 

Given a graph G — (V, E) and a Hamiltonian H acting 
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FIG. 2: Left: Local Hamiltonian on a 5-vertex graph. Each 
edge in the graph is a clique, but the associated Hamilto- 
nian terms hg — are null. The four triangular cells in 
the graph are three-body cliques, and the associated Hamil- 
tonian terms are h v = o\a\a\, — a^cr^a^, his = cr%o%<j\, 
ho = a%G\(J%. Right: Lattice tiling with the graph on the 
left. 

on the particles located in V , we say that H is shield 
commuting when it can be written as H = Hab + Hbc 
with [Hab, Hbc] = whenever B shields A from C . 

Thus, Eq. Q implies a complete equivalence between 
positive quantum Markov networks and Gibbs states aris- 
ing from shield commuting Hamiltonians. Clearly, the lo- 
cal commuting Hamiltonians form a subset of shield com- 
muting Hamiltonians. In particular, it can be easily ver- 
ified that the terms of a local Hamiltonian commute, but 
the shield commutation condition involves many-body 
operators, and there is no obvious efficient way to test 
it. 

It is the purpose of the this section to examine whether 
shield commuting Hamiltonians are equivalent to local 
commuting Hamiltonians. In ID it follows from the 
decomposition given in Eq. [2] that a positive Markov 
state \ogp = H may be decomposed as H = J^hii+i. 
Equation Q then implies that for each pair of terms 
[hi-u, ^ii+iT — 0- Such an equivalence has also been 
shown for cycle- free graphs in |23j , and in Theorem [4] we 
will extend it to all graphs containing no triangles. For 
more general lattices however, it is not always the case 
that a shield commuting Hamiltonian is a local commut- 
ing Hamiltonian, as will be demonstrated by the following 
example. 

We now establish that there exist positive quantum 
Markov networks that are not Gibbs states of local com- 
muting Hamiltonians. Consider the graph G illustrated 
on the left of Fig. [2] and its associated Hamiltonian 
H = hy + h < + ft, a + ^> • There are only two choices of 
regions A, B, and C for this graph such that B shields 
A from C, namely A = {1}, B = {2, 4, 5}, C = {3}; and 
A = {2}, B = {1,3,5}, C = {4}. For the first choice, we 
have Hab = h v + h > and Hbc — ' l <i + ^A > and we can 
easily verify that [Hab, Hbc] — 0, which ensures that 
I (A : C\B) = 0. The second choice has Hab = h\7 + ^< 
and Hbc = h> + h&, an d leads to the same conclu- 
sion. Thus, the pair (p = ^e H , G) is a positive quantum 
Markov network. Nevertheless, p is not the Gibbs state 



of a local commuting Hamiltonian as the individual terms 
in the Hamiltonian do not mutually commute. 

The choice of the Hamiltonian leading to this situation 
is very contrived, and the presence of triangular cliques 
is crucial in our construction. In fact, we can rule out the 
existence of such examples in the absence of triangular 
cliques. 

Theorem 4 Let G = (V, E) be a graph and p be a density 
matrix for the particles located at the vertices of G. If the 
pair (p, G) is a positive quantum Markov network and G 
contains only two-body cliques, then the state p can be 
expressed as p = e where 

H=J2h Q , [h Q ,h Q ,} = (10) 
Qec 

is the sum of mutually commuting Hermitian operators 
hq on the particles located in cliques Q. 

This theorem can be seen as a generalization of a The- 
orem presented in the supplementary material of (23) . 
which is restricted to cycle-free graphs. 

Theorem 5 (Poulin and Hastings |23j) Let G — 

(V, E) be a tree graph and p be a density matrix for the 
particles located at the vertices of G. The pair (p, G) is a 
positive quantum Markov network if and only if the state 
p can be expressed as p = e H where 

H=J2 h Q> [h Q ,h Q ,]=0 (11) 
Qec 

is the sum of mutually commuting Hermitian operators 
hq. on the particles located in cliques Q. 

The "if" implication is a corollary of Theorem [3j To 
prove the "only if" direction of Theorem |4j we first need 
to introduce the concept of a genuine fc-body operator, 
closely related to the cumulant expansion. Let Hx be an 
operator on some set of particles X = {1,2, ... ,k}. We 
say that Hx is a genuine operator on X if TiyH = 
for all non-empty Y C X . This implies that the cumu- 
lant expansion of Hx contains only one term, Hx itself. 
Equivalently, we can define a cumulant expansion as the 
decomposition of an operator into a sum of genuine op- 
erators on all subsets of particles. One useful fact about 
genuine operators is given by the following lemma. 

Lemma 1 Let Hab be a genuine operator on AB and 
Hbc be a genuine operator on BC. Then [Hab, Hbc] * s 
either or a genuine operator on AB'C where B' C B 
and B' £ 0. 

This lemma is proven by considering the operator- 
Schmidt decomposition of the two operators 

H AB =Y,F j A ®G j B (12) 

3 

H bc = Y. R b® S c (13) 

k 
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where {F^}, {G B }, {R B }, and {S c } are sets of orthog- 
onal operators. Then, 

[Hab,H bc ] =Y. F a® \G^R k B ] ® Sc. (14) 

jk 

The lemma follows from the fact that the commutator 
[G g) Rg], if non-zero, cannot be proportional to the iden- 
tity for finite dimensions. 

We now turn to the proof of Theorem |4| In the case 
of a graph with only two-body cliques, Theorem [2] al- 
ready establishes that the cumulant expansion of the 
Hamiltonian is a sum of terms on edges and vertices, 
H = E ee£ K e + ^2 vev K v . We will prove Theorem[4]by 
expressing the Hamiltonian as H — J2 e eE ^ e + ^2vev ^ v 
where all terms mutually commute. Note that the terms 
h e and h v can differ from the cumulants K e and K v . In 
particular, a one-body cumulant K v can be split arbi- 
trarily among h v and the edge terms h e having vertex v 
as an endpoint, so h e needs not be a genuine operator on 
the edge e. 

We first prove that all two-body cumulants K e mu- 
tually commute. Pick any disjoint regions A, B, and 
C = V — A — B, such that B shields A from G, so we 
can write H = Bab + H BC with [Bab, Hbc] = 0, c.f. 
Eqs. (|3j [2]). We can express Bab and Bbc in terms of 
the cumulants of B, with some cumulants contributing 
to Bab and some to Bbc- The contribution of most 
cumulants is unambiguous, but the ones supported only 
on region B could be either attributed to Bab or Bbc- 
More generally, the cumulant expansions have the form 

B AB =J2 Kab + Y. Ka + l>2 K ™' +J2 K x ( 15 ) 

ab a aa' XcB 

h bc ^Y, Kb - + Y, K - + Yl K ^ + S K * (16) 

be c cc' X(ZB 

where 

K*+K c x = K w +Y,K h . (17) 

XcB bb'eB b 

and where it is understood that small case letters run 
only over the vertices in the set labeled with the same 
capital letter, i.e. a g A, b € B, and c G C. In these 
equations, we use a double index such as ab to denote 
an edge e = (a, b) and set K a b = if (a, 6) ^ E. In 
those terms, the commutation relation [Bab, Hbc] — 
becomes 

a be 

bc,XcB ab,XcB 
X,YCB 

Note that each summand in the first line has a support 
that differs from all other terms in the equation. By 




FIG. 3: Definition of the regions A, B, and C for a given 
neighbor v\ of vertex u. 

Lemma[lj they must be individually equal to 0. By vary- 
ing over all possible choices of regions A, B, and C, the 
first line shows that the two-body cumulants mutually 
commute, as claimed. 

Now that we have established that the two-body cu- 
mulants K e mutually commute, our goal is to split each 
single-body cumulant K u among the edge terms h e act- 
ing on site u and the vertex term h u to obtain mutually 
commuting terms. More formally, let Af(u) = {v € V : 
(u, v) 6 E} denote the neighborhood of u, and denote 
the degree of u by d(u) = |A/"(u)|. We will decompose 
K u = h u + J2veAf(u) G«) and define the edge terms 

h U w = Kuw + G w + (18) 

in a way that [h u , h uv ] = and [h uw , h uv ] = 0, or equiv- 
alently 

[h U w, h uv ] — [K uw + G w + G^, K uv + G" + G u ] (19) 

= [K U w, G v u ) + [G™, K uv ) + [G u , G™] = 0. 

(20) 

Since the three commutators in the last line have distinct 
support, by Lemma [l] this last condition is equivalent to 
[Kuw,G v u ] = and [G V U ,G^] = for all (u,w) € E and 
(u, v) E E. What remains to be explained is how to 
choose the G u . 

Consider a vertex u in the graph. Let {v\, V2, ■ ■ -Wfe} = 
J\f(u) denote the immediate neighbors of u. Define three 
sets of vertices 

AcN{u) (21) 
B = AT (A) (22) 
G = V - A - B. (23) 

where the neighbors of the set A are defined in the natural 
way N(A) = {w G V : 3v G A,(w,v) G E}. These 
definitions are illustrated on Fig. [3] In other words, B 
consists of the neighbors of a subset A of the neighbors 
of u, so in particular B includes the vertex u itself. It 
follows that B shields A from G. Moreover, because the 
graph has only two-body cliques, all the neighbors vi ^ A 
are contained in G. The Markov condition implies that 
H = Bab + H B c with [B AB , B B c] — for this choice 
of regions, c.f. Eqs. ([3j[4|. 
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Consider the cumulant expansion of Hab and Hbc- 
Because H = Hab + Hbc-, the cumulants in these ex- 
pansions can only differ from those of H on the region 
B, i.e. we must have 



H AB =J2 K *b + J2 Ka + Yl K * ( 24 ) 

ab a XCB 

Hbc = Kbc + J2 K - + £ K x + J2 K ™' ( 25 ) 



be 



XCB cc' 



where 



£ K$+ Kx = £ K bb>+J2 Kb - ( 26 ) 



XCB XCB bb'CB-u 



Note that Eq. (24) cannot contain two-body cumulants 
K aa ' because this would create a triangle in the graph, 
which is ruled out in the hypothesis of the Theorem. For 
the same reason, the two-body cu mula nts Kw in Eq. ( 26 ) 
cannot include node u. Equation 26 implies that K = 
—Kx for all \X\ > 2 as well as for all \X\ = 2 when 
u e X. Now, consider the commutation [Hab, Hbc] = 
in this expansion: 



= J2[Kab,K bc ] 

abc 

+ Y / [Kab,K^]+Y / [Kx,K, 



abX 



bcX 



(27) 
(28) 
(29) 



XY 



The first line Eq. ( 27 1 is since two-body cumulants of H 
mutually commute. By Lemma [l] each individual term 
in the second line must be because they each have a 
distinct support from all other non-zero terms in the sum. 
In the last line, by Lemma [T] the only terms in the sum 
that can be supported on the single site u are of the 
form [K£,K°] or [K$,K%] with u e X and \X\ > 1, 
but these latter are since = —K^ as shown above. 
Therefore, the terms with X = Y = u are decoupled from 
the other terms, and so must obey [K^,K^] = 0. We 
conclude that for any subset A C W(u), we can express 
K„ = K A + K?. with 



[K£,K%] = [K au ,KS] = [KtK u 







(30) 



for all a € A, and c e Af(u) - A. 

Consider the subgraph G{u) of G consisting of site u 
and its neighbors Vj, and define the Hamiltonian H(u) = 
~Y1, V K UV + K u . Note that G(u) is a tree and that the 
pair (p(u) = |e fl W,G(ti)) is a positive quantum Markov 
network. This last fact follows from the commutation 
relations of Eq. ( 30 ) and Theorem [3j By Theorem [5] 
it follows that K u can be decomposed into K u — h u + 
J2veM(u) G u in sucri a way that h' uv = K uv + G v u and 
h u all mutually commute. Repeating the argument for 



every node u € V and defining h v 
complete the proof of Theorem [4] 



K -X- C v 



Gu 



Note that the proof of Theorem [4] makes crucial use 
of the fact that the lattice has only two-body cliques. 
As a consequence, the proof does not extend directly to 
Hamiltonians that are the sums of two-body terms but 
embedded on, say, a triangular lattice. 



IV. CONCLUSION AND DISCUSSION 

In this article, we have presented a simplified proof 
that all positive quantum Markov networks living on lat- 
tices of finite dimensional quantum systems are Gibbs 
states of Hamiltonians local to the cliques of the cor- 
responding graph; proved that there is an equivalence 
between Gibbs states of commuting Hamiltonians and 
positive quantum Markov networks on graphs that do 
not contains triangles, shown that Gibbs states arising 
from local Hamiltonians with commuting terms are posi- 
tive quantum Markov networks; and lastly; demonstrated 
that there exist positive quantum Markov networks on 
triangular lattices that cannot be written as Gibbs states 
of commuting local Hamiltonians. While both classical 
and quantum positive Markov networks are Gibbs states 
of Hamiltonians local to the cliques of their underlying 
graph, the essential difference is that the terms of a clas- 
sical Hamiltonian commute by construction, whereas the 
terms in the quantum Hamiltonians are only required to 
commute when grouped into globally defined regions as 
in FigfTJ 

It remains an open question whether there exist Hamil- 
tonians that satisfy this global commutation property — 
that we called shield commutation — but that cannot be 
transformed into a local commuting Hamiltonian under 
a suitable renormalization procedure. Despite many at- 
tempts, we have been unable to construct such models. 
This is illustrated by the example of Sec. |III C| If the 
network is extended to a regular 2D lattice as shown on 
the left of Fig. [2j we again obtain a Gibbs state that is a 
quantum Markov network, yet the Hamiltonian terms do 
not commute. However, if we coarse grain the lattice by 
combining the central spin of each unit cell to the spin 
immediately to its northeast — and thus join /i v and 
and similarly h& and to the same cliques — we obtain 
a local commuting Hamiltonian. This illustrates that the 
lattice model can be fixed by a local rearrangement of the 
degrees of freedom, but has no "large-scale" obstructions 
to commutation. 

It has been established [H H] that for Hamiltonians 
containing only two-body commuting interactions, the 
Hilbert space of each vertex must split into a direct sum 
of factor spaces as in Eq. [2] This is reminiscent to the 
fact that we can only prove an equivalence between local 
commuting Hamiltonians and positive quantum Markov 
networks in the absence of triangular cliques. Although 
the local splitting property is not necessary for a fac- 
torization into commuting terms — the projector into the 
code space of Kiteav's toric code [12] is a quantum 
Markov network which factorizes by construction — , it is 



sufficient. Since commutation of genuine operators that 
act non-trivialiy on more than two subsystems does not 
imply that the local Hilbert spaces will split [3] , quantum 
Markov networks containing triangles and hence three- 
body cliques may involve complex non-local structures, 
which could allow for the shield and local commuting 
properties to be inequivalent even under coarse graining. 

A folk theorem in condensed matter physics states that 
all phases of matter can be realized with local com- 
muting Hamiltonians. This is supported by the fact 
that Markov networks are fixed points of a renormal- 
ization procedure. A corollary of this statement would 
be that, under a suitable renormalization procedure, all 
Gibbs states are quantum Markov networks, thus es- 
tablishing a complete equivalence between Gibbs states 
and quantum Markov networks. If models with large- 
scale obstructions to commutation could be found on 
the other hand — Hamiltonians that are shield commuting 
but cannot be locally transformed into local commuting 
Hamiltonians — , they would reveal a new phase of matter 
that exhibits quantum non-locality without long range 
entanglement [3] , and would require a refinement of this 
folk theorem. They would also form a class of Hamiltoni- 



ans of intermediate complexity between commuting and 
general local Hamiltonians pj IH [TTJ [24] . 

Lastly, the relation between the Markov condition 
and area laws is also intriguing. At finite temperature, 
the entropy of a region should be mostly extensive — 
scaling with its volume — but the zero-temperature area- 
law should translate in an additional boundary contri- 
bution. Thus, we can expect a generic scaling S(A) = 
a\A\ + f3\dA\ for large enough regions A. A simple ge- 
ometric argument shows that under this scaling, regions 
A, B, C chosen as in Fig.[l]always obey the Markov con- 
dition. This gives us additional reasons to believe that 
Markov networks are generic properties of Gibbs states 
on sufficiently coarse grained networks. 
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